Anisotropic flow in striped superhydrophobic channels 
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We report results of dissipative particle dynamics simulations and develop a semi-analytical theory 
of an anisotropic flow in a parallel-plate channel with two superhydrophobic striped walls. Our 
approach is valid for any local slip at the gas sectors and an arbitrary distance between the plates, 
ranging from a thick to a thin channel. It allows us to optimize area fractions, slip lengths, channel 
thickness and texture orientation to maximize a transverse flow. Our results may be useful for 
extracting effective slip tensors from global measurements, such as the permeability of a channel, in 
experiments or simulations, and may also find applications in passive microfluidic mixing. 



I. INTRODUCTION 

Textured surfaces play a major role in microfluidics, 
since the high surface to volume ratio enhances fluid- 
surface interactions An important class of phe- 
nomena involve "transverse" hydrodynamic couplings in 
anisotropic channels, where an applied pressure gradi- 
ent or shear rate in one direction generates flow in a 
different direction, with a nonzero perpendicular com- 
ponent. Transverse hydrodynamic couplings in pressure- 
driven flow through a textured microchannel were ana- 
lyzed theoretically [l-Ill and applied to passive chaotic 
mixing in a herringbone grooved channel 0, Q . Such mi- 
crofluidic devices have also recently been used to separate 
or concentrate suspended particles 0. 

Here we deal with a special (and different) type of sur- 
face textures, namely a superhydrophobic (SH) surface in 
the Cassie state, where air micro- and nanobubbles are 
favored, and can generate a number of amazirig proper- 
ties [1] , including a very large liquid slippage d, hfl. Such 
a slip also occurs at smooth hydrophobic surfaces jl ll [l^ . 
but with a relatively low amplitude, characterized by a 
slip length (extrapolated distance on which the liquid ve- 
locity vanishes) of the orders of tens of nm or smaller 
[TgI . Due to a local slip length of tens of ^m over gas re- 
gions the dramatic enhancement of slip at such a super- 
hydrophobic surface could be achieved [T7l - l20| and may 
dramatically reduce viscous drag in a channels. 

Directional SH textures generate the anisotropy of ef- 
fective slip, which becomes a tensor, boff = {bf^}, repre- 
sented by a symmetric, positive definite 2x2 matrix 




FIG. 1: Sketch of the symmetric striped channel (a): O = 
7r/2 corresponds to transverse stripes, O = to longitudinal 
stripes; (b) situation in (a) is approximated by a periodic cell 
of size L, with equivalent flow boundary conditions on the 
gas-liquid and solid-liquid interfaces. 



a given geometry. For all anisotropic surfaces its eigen- 
values and hj:^^ correspond to the fastest (greatest 
forward slip) and slowest (least forward slip) orthogonal 
directions [Sl- Note that the concept of an effective slip 
length tensor is applied for an arbitrary channel thick- 
ness 21 1 , being a global characteristic of a channel [l^l ■ 



diagonalized by a rotation with angle Q 

„ ( cos 6 sin O 
~ I -sine cose 



(1) 



(2) 



The tensor formalism allows us to easily change the ori- 
entation of a texture, once a problem has been solved for 



A corollary of this is that the eigenvalues depend not only 
on the parameters of the heterogeneous surfaces (such as 
local slip lengths, fractions of phases, and a texture pe- 
riod), but also on the channel thickness. However for a 
thick (compared to a texture correlation length) chan- 
nel they become a characteristics of a heterogeneous in- 
terface solely 0. In the general case of any direction 
e, the anisotropy of effective slip means that the flow 
past such surfaces becomes misaligned with the driving 
force. Therefore, anisotropic SH textures can potentially 
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be used to generate transverse hydrodynamic flow, which 
is of obvious fundamental and practical interest. How- 
ever, most of the prior theoretical work has focused on 
the maximization of the effective slip, and general prin- 
ciples to optimize transverse phenomena have received 
much less attention and have not yet been established. 

Computer simulations might be expected to shed some 
lights on a flow past anisotropic surfaces. Earlier Molec- 
ular Dynamics (MD) simulations observed how modula- 
tions of hydrophobicity induce (small) variations in local 
slip at channels of finite thickness, but did not attempt to 
relate these to effective boundary conditions [l^, [23|. In 
efforts to better understand the connection between the 
parameters of the texture and the effective slip lengths, 
several other groups have been performed MD simula- 
tions of the fiow past striped anisotropic surfaces. Most 
of these studies considered a thick channel geometry, and 
focused only on calculating eigenvalues of the effective 
slip- length tensor [l^, HH • There are however two sim- 
ulations which are directly relevant. Recent work (26j 
studied the dependence of a ratio of transverse and for- 
ward effective slip velocities on the flow orientation in 
a thick channel, but it did not explore the situation of 
a finite channel thickness. It has been found that the 
theoretical predictions based on a concept of tcnsorial 
slip consistently overestimated the results of MD sim- 
ulations, which has been attributed to the fact that the 
stripe width was comparable to the molecular size. Other 
recent simulations [2l[ used the lattice-Boltzmann (LB) 
approach to calculate a fiow in an asymmetric channel 
(i.e. with one SH and one no-slip hydrophilic wall) of ar- 
bitrary thickness. This study confirmed that the eigen- 
values of the effective slip-length tensor depend on the 
gap and validated the concept of tensorial slip. An ex- 
cellent agreement of LB results with a continuous theory 
gave a strong support to macroscopic arguments. How- 
ever, no attempts have been made to simulate the de- 
fiection of fiow from the main direction, which for an 
asymmetric SH channel was not expected to be very dra- 
matic [Toj . 

The goal of the present paper is to provide some gen- 
eral theoretical results to guide the optimization of trans- 
verse hydrodynamic phenomena in a symmetric channel 
with two aligned striped SH walls. In the case of thin 
channels, symmetric striped walls provide rigorous up- 
per and lower bounds on the effective slip over all pos- 
sible two-phase patterns [l^l and are expected to gener- 
ate very strong transverse fiow [28[, which could lead to 



very efficient mixing upon spatial modulation of the tex- 
ture. Here we consider a channel of arbitrary thickness, 
and optimize its orientation to maximize pressure-driven 
transverse fiow. Our consideration is based on theoretical 
analysis and dissipative particle dynamics (DPD) simu- 
lations. By using relatively wide stripes as compared 
with a recent MD study [l^ we find that our results in a 
very good quantitative agreement with the macroscopic 
theory. 



II. MODEL AND GENERAL CONSIDERATION 

We consider a channel consisting of two equal symmet- 
ric SH walls located at y = and y = H and unbounded 
in the x and z directions as sketched in Fig. [1] The origin 
of the coordinates O is placed in the plane of a liquid-gas 
interface at the center of the gas sector. The x axis is 
defined along the pressure gradient. As in previous pub- 
lications d, [H, m, [2^ HO] , we model SH plates as fiat 
interfaces with no meniscus curvature, so that they ap- 
pear as being perfectly smooth with a pattern of bound- 
ary conditions. The latter are taken as no-slip (61 ~ 0) 
over solid/liquid areas and as partial slip (62 ~ b) over 
gas/liquid regions. In this idealization, by assuming a fiat 
interface, we neglect an additional mechanism for a dis- 
sipation connected with the meniscus curvature [sil [3^ . 
The fraction of the solid/liquid areas is denoted as (/>!, 
and that of the gas/liquid areas as (/)2 = 1 — ^1. 

The flow is governed by the Stokes equations 



Vp, V • V = 0, 



(3) 



where u is the velocity vector, and the average pressure 
gradient is always aligned with the x-axis direction: 



(Vp) = (-a,0,0) 



(4) 



The local slip boundary conditions at the walls have the 
form 

v(x,0,z) = 6(a;,z) • — (a;,0,z), y • v(a;, 0, z) = 0. (5) 

In addition, a symmetry condition applies at the mid- 
channel. 



9v 

dy 



{x,H/2,z) = 0, y.v(.T,i7/2,z)=0. (6) 



Here the local slip length b{x, z) at the SH surface is 
generally a function of both lateral coordinates. 

The effective slip length b^s at the SH surface is deflned 
as usual. 



fccff 



(7) 



where (. . .) denotes the average value in the plane xOz. 
It can equivalently be defined via a permeability tensor 



(8) 



In the linear response regime, the averaged flow rate, 
(Q), is proportional to (Vp) via the permeability tensor. 



(Q) = — k.(Vp), 



(9) 
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which may be rewritten as 

(Q>x = ^ (fcllcos^e + fc-Lsin^e) , (10) 

(g)^ = - (^fcii - /j-l) sine cose. (ii) 

Integrating the velocity profile across the channel we 
obtain 



H 



(Q) 



(v(y)) dy, 



H 



{vx{v))dy, 



H 



{vz{y))dy. 



(12) 

(13) 
(14) 



The above formula are valid regardless of the thickness 
of the symmetric SH channel and arc independent of the 
details of the textured surface. There could be arbitrary 
patterns of local slip lengths, and the latter could itself 
be a spatially varying tensor. 



III. THEORY FOR STRIPED PATTERNS 

To illustrate the general model, in this section we focus 
on flat patterned SH surfaces consisting of aligned peri- 
odic stripes, where the local (scalar) slip length b varies 
only in one direction. We mostly follow the approach 
developed before [2l|, but apply the method to a sym- 
metric channel situation as sketched in Fig. [1] where the 
surfaces are covered with arrays of gas/liquid stripes with 
width S and period L. Thus, (f>i = l — S/L and 02 = ^/L, 
respectively. 

For transverse stripes, the flow is two-dimensional 
V = {v^X^,y),Vy{x,y),0), V:,{x,0) = b{x)dyVx{x,0), and 
fj,(x, 0) ~ 0. For longitudinal stripes, we also have a 
plane flow: v = {vx{y, z),Q,Q), Vx{Q,z) = b{z)dyVx{0, z). 
As the problem is linear in v, we seek for a solution of 
the form 



.(1) 



(15) 



where v^^^ is the velocity of the usual no-slip parabolic 
Poiseuille flow 

^^"^=[^y(H~y),o,oy (16) 

and v^^) is the SH slip-driven superimposed flow. 



A. Longitudinal stripes 

In this situation the problem is homogeneous in x- 
direction {d/dx = 0). The slip length b{x,z) = b{z) 



is periodic in z with period L. The elementary cell is de- 
termined as b{z) = & at \z\ < S/2, and b{z) = at S/2 < 

\z\ < L/2. In this case the velocity v^^'' = {vi^\ 0, 0) has 
only one nonzero component, which can be determined 
by solving the Laplace equation 



w\i'Hy,z) = Q, 



(17) 



with the conditions ((S])-®. 

The Fourier expansion of a periodic solution satisfying 
([6]) reads 



vi'^ (y, ^) = y + E cos(A„z)e-^"^[l + e'-^^v-"^] 

71=1 

(18) 

with A„ = 21m/ L. The sine terms vanish due to symme- 
try. Applying ([5]) we then obtain a trigonometric dual 
series: 



^ + E a„ [1 + b\n tanh(A„i//2)] cos(A„z) = 6, (19) 

\z\ < 6/2, 



n=l 



^ + ^a„cos(A„z) = 0, 5/2<\z\<L/2, (20) 



where 



ao 



2?/ao 



2r/a„(l 



(tH 



n> 1. 



The dual series , (pOj) provide a complete description 
of the hydrodynamic flow and the effective slip length 
big = olq/2 (due to Eq. ([T2|)) in the longitudinal direc- 
tion, given all the stated assumptions. These equations 
can be solved numerically, but exact results can be ob- 
tained in the limits of thin and thick channels. 

For a thin channel, <C -L, we can use that 
tanh<|j_j_g = 0(t). By substituting this expression into 
(|19p and keeping only values of the first non-vanishing 
order, we find 



5/2 
-5/2 



(21) 



This is an exact solution, representing a rigorous upper 
Wiener bound on the effective slip over all possible two- 
phase patterns in a thin channel 27j . 

In the limit of a thick channel, 7? 3> i, we can use 
that tanh(i — > oo) — >■ 1 and the dual series (fT9 |) - (f20| can 
be solved exactly, giving [s^ 



-'efl 
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+ tan 







(22) 
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B. Transverse stripes 

To describe a two-dimensional flow in this case we use 
a standard technique and introduce a stream function 
ip{x, y) and the vorticity vector u)(x, y). Thus, the veloc- 
ity field is represented by v(a:, y) = {dtjj/dy, —dtp/dx, 0), 
and the vorticity vector, u:{x, z) = V x v = (0, 0, w), has 
only one nonzero component, which is equal to 



(23) 



The solution can then be presented as the sum of the base 
flow with homogeneous no-slip condition and its pertur- 
bation caused by the presence of stripes as 



(7 y'^ (tH 
^ 7r + ^ + 

a (jH 
w = -y - — + a;i. 

The problem for -i/)! and loi reads 

vVi = -wi, v^wi = 0, 



(24) 
(25) 

(26) 



which can be solved with respect to the BC ([S])-® and 
an additional condition that reflects our definition of the 
stream function: 



^-1(2^,0) = 
The general solution reads 

i^i{x,y) = Poy 



(27) 



00 



00 



/ J \ n 




e "^cosA„x 



2 



^ a„ cos(A„x) = 0, S/2<x<L/2, (31) 



n=l 



Here, 



2 



2ri 



sinh(A„-ff) 



and 



Wit) 



H 



cosh(t) — 1 



(32) 

P„, n > 1, (33) 



sinh(t) — t 

In the limit of a thin channel one has W{t)\t^ao 
Zt^^ + 0{t), from which it follows that 



(34) 



bH4)2 



(35) 



H + Qh(t)i 

This exact equation represents a rigorous lower Wiener 
bound on the effective slip over all possible two-phase 
patterns in a thin symmetric channel 27 1 . 

For completeness we discuss again the two limiting sit- 
uations: 

Ocs\H<^b.L 



6. 



■-H oc H, 



OC 6. 



(36) 
(37) 



In the thick channel limit, the dual series (|30l) and 
(PT|) take the same form as in prior work [l^ (due to 
W{x — )- 00) — )- 1), whence we derive [11] 



^ctf 



L 




In 








-Am 

27rfe 


sec 4- tan 



(38) 



UJi 



00 

{x,y) = J2 + Mi^^e-^-y) cos(A„x). 



(29) 

The following relations between Fourier coefhcients may 
be established: 

p(i) - _P(2) = _p 



M^'^ = -^sinh(A„i//2)e-^"^/2p„, 



M(2) = l^sinh(A„i7/2)e^"^/2p„. 

Again we obtain a dual series problem, which is similar 
to (HI]) and (HOI): 



^ + ^ a„ [1 + 26A„W^(A„iI)] cos(A„x) = b, (30) 

0<x< 6/2, 



C. Transverse flow 

The tensorial nature of the effective slip is physically 
due to secondary flows transverse to the direction of the 
applied pressure gradient. Now we focus on transverse 
flow optimization, which is necessary for a passive mixing 
in a symmetric SH channel. Our aim is to optimize the 
texture, channel thickness and the angle & between the 
directions of stripes and the pressure gradient, so that 
\{Qz) / {Qx)\ is as large as pos sible. 

Following the approach [lOl , one can derive 



F = 



6^) sin 9 cos 



H + 6blff cos2 e + 66^ff sin^ 8 



(39) 



Consider now the tilt angle, Q, confined between and 
tt/2 and rewrite Eq. ([5^ as 



6ie - l)r 



(40) 
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where r = tan 9, = bl^/h^^ and /i^ = H/h^^^. By 
evaluating dF/dr = 0, we find that the maximum occurs 
at 



T, — tan 0* = 



1/2 



and its value is 

F(r,) 



He - 1) 



[(/l2+6e2)(/l2+6)]l/2 



(41) 



(42) 



In the limit of a thick channel, H ^ L, h ^ oo (owing 
to the fact that 6^ is limited by the local slip length h 
and independent of H for that case), and Eq. (|4T|) gives 



tan0»|^^^ — 1; 



t/4, 



(43) 



and, correspondingly. 



He 1) 



H/L- 



= 0{L/H). (44) 



Therefore, the mixing in a thick symmetric SH channel 
would be not very efficient. According to our formula, 
to maximize \{Qz) / {Qx)\ an efhcient strategy would be 
to use striped textures with largest physically possible b 
and a very low fraction of a solid phase, 02 — >■ 1. 

When the channel is thin, H L, the interplay be- 
tween b and H gives two possibilities. First, if & <C -ff ^ 
L, the slip length is isotropic, &^ = 602(1 — Q(t>ib/H) ~ 
and /i^ = H/{b(j)2) > 1, ^ ~ 1. Therefore, again 
0, ~ 7r/4, yet the value 



F(r,) 



i72 



= O 



n 2^ 



(45) 



is negligible. Then, for H <C min{6, L} we get 
~ 602/01 = const, 

= 660i/H » 1. 

Substitution of these expressions into Eqs. (|4T|) and (|42 
gives 



tan0 



F(t,) 



*lif<Smin{6,L} 



1 



H 



1/2 



(46) 



Lff<min{b,L} ~ 2 \ H 



1/2 



> 1. (47) 



Since 0i02 = (1 — 02)02 is maximal for 02 = 0i — 0.5, 
the direction of optimal inflow angle is 

which coincides with the results of [l^ derived by using 
a different method. 



IV. SIMULATION METHOD 

Our simulations are done using Dissipative Particle 
Dynamics (DPD), an established method for mesoscale 
fluid simulations, which is fully off-lattice and particle 
based and naturally includes thermal fluctuations [33 - 
l36j . More specifically, we use a DPD version without con- 
servative interactions, and combine that with a tunable- 
slip method that allows one to implement arbitrary hy- 
drodynamic boundary conditions |37|. The detailed im- 
plementation of the tunable-slip method can be found in 
Ref. [13 ■ In the following we only give a brief description 
and introduce the simulation parameters. 

In the tunable-slip boundary approach, the interaction 
between the channel walls and the fluid particles has two 
contributions. The first is a Weeks-Chandler- Andersen 
(WCA) interaction to mimic the impermeable surfaces. 



y < 2i/V 

y > 2i/V 



(49) 



where y is the distance between the fluid particle and 
the wall. This is just a Lennard- Jones interaction with 
a cutoff at the potential minimum, corresponding to the 
pure repulsive part of the potential. In the following, 
the WCA parameters will set our simulation units, i.e., 
the energy unit e and the length unit a. The third unit 
is the particle mass m. The second part of the wall- 
fluid interaction is a coarse-grained friction force, which 
is introduced in a similar spirit than the DPD approach. 
Specifically, the effect of the wall friction on i-th particle 
is implemented by introducing a pair of Langevin-type 
forces, 



The dissipative contribution has the form 
Ff = -7La;L(2;)(vj - v^^aii)- 



(50) 



(51) 



This force is proportional to the relative fluid velocity 
with respect of the wall, and the proportionality factor 
is the local viscosity ^llul(z). The position-dependent 
function ujl{z) is a monotonically decreasing function of 
the wall-particle separation, and vanishes when the fluid 
particle is further away from the wall than a cutoff dis- 
tance Zc. In our simulations, we take u>l{z) to decrease 
linearly from 1 to 0. The prefactor 7^ characterizes the 
strength of the wall friction and can be used to tune the 
value of the slip length. A random force obeying the 
fluctuation-dissipation theorem is required to ensure the 
correct equilibrium statistics. 



(52) 



where each component of is a Gaussian distributed 
random variable with zero mean and unit variance. 

The simulations are carried out using the open source 
simulation package ESPResSo [l^. All simulations are 
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performed with a time step At = O.Ol^m/ecr, and the 
temperatm'e of the system is set at ksT = le. The fluid 
has a density p = 3.75a~^. Fluid particles have no con- 
servative interactions, they interact only with the dissipa- 
tive part of the DPD interactions, The DPD interaction 
parameter is chosen at jdpd ~ 5.Qy/me/a and the cut- 
off radius is l.Ocr, which results in a shear viscosity of 
r]s = 1.35±0.01v^/cr2. 

The slip length b can be calculated analytically as a 
function of the simulation parameters (the wall friction 
parameter and cutoff Tc) to a very good approxima- 
tion. For the purpose of this work, however, we need 
highly accurate values for both b and the position of the 
hydrodynamic boundary, therefore we determined it with 
the method described in Ref. [37| by simulating plane 
Poiseuille flow and plane Couette flow. Fig. [5] shows the 
relation between the slip length and the wall friction pa- 
rameter 7i, for a wall interaction cutoff Vc = 2.0a. By 
setting the wall friction parameter 7^, we can adjust the 
slip length to arbitrary values. The no-slip boundary 
condition is implemented with 7^ = 5.26y/rne/a. The 
position of the hydrodynamic boundary, was determined 
to be 1.06 ± 0.12(7 away from the simulation wall. 



350 




0.01 0.1 1 10 



FIG. 2: The relation between the slip length 6 and the wall 
friction parameter jl, for a wall interaction cutoff 2.0a. The 
inset shows an enlarged portion of the region where the slip 
length is zero. The no slip boundary condition can be im- 
plemented by using = 5.26^me/a. Dashed curves show 
analytical predictions [stI ]. 

For simulations of patterned surfaces, we use a box 
of 20a X (H + 2a) x 50cr. The two impermeable sur- 
faces lie parallel to the xz-plane and are separated by 
a distance iJ -I- 2a. Periodic boundary conditions are 
assumed in the xz plane in the directions parallel and 
perpendicular to the stripes. We have tested that in- 
creasing the simulation box does not change the out- 
come. The surfaces are decorated by alternating no-slip 
and partial slip stripes. The pattern has a periodicity 
of L = 50(7. The wall friction parameters are chosen 
7l = 5.26\/?7ie/(7, which implements a no-slip boundary 
condition, and 7^ = 0.02820-y/me/(7, corresponding to 



b = 50(7. A simulation system is completely specified by 
the following parameters: (L, 6, (/)2, H^- 

The simulation starts with randomly distributed DPD 
particles. An external body force is assigned to each par- 
ticle to mimic the pressure gradient. The amplitude of 
the force is adjusted such th at the m aximum shear rate 
at the surface is 7 = O.Oli/m/ecr^. Small shear rates 
are necessary to avoid shear-rate dependencies of the slip 
length [3^ and reduce inertia effects. With these param- 
eters, typical Reynolds numbers in our system are still of 
order 10, which is much larger than typical values in mi- 
crofluidic setups. To reach realistic Reynolds numbers, 
one would have to reduce the body force by four orders 
of magnitude. Unfortunately, the necessary simulation 
time for gathering data with sufficiently good statistics 
would then increase prohibitively, since the values of the 
relevant observables become very small and the correla- 
tion times increase with 1/7. Test runs with selective 
parameters were performed using a smaller body force, 
but the results did not change significantly. Therefore, we 
chose to work at these relatively high Reynolds numbers 
(we are still deeply in the laminar regime), and tolerate 
slight inertia effects, which we will discuss further below. 

One important criterion for a correct simulation of the 
pressure-driven flow is a constant density profile. Fig. 
[3ja) shows an example of the fluid density profile for 
p = 3.75(7^"^. The channel has a thickness of i7 = 50(7. 
The periodicity of the striped pattern is L = 50(7, and 
the partial-slip part has a slip length oib ~ 50(7 and area 
fraction (\)2 — 0.5. The fiuid density is uniform in the 
middle of the channel, and drops to zero in the prox- 
imity of the wall due to the repulsive interaction of the 
wall on the fluid particles. This leads to a physical wall 
position which is approximately 1.0(7 away from the sim- 
ulation wall. The quoted value of channel thickness H is 
the separation between the physical walls. In contrast to 
earlier MD simulations [25| , there are no molecular layer- 
ing effects and density oscillations near the walls, because 
our DPD particles have no conservative interactions. 

Fig. [2Jb) shows the averaged velocity profile for the 
same system. The system is allowed to run up to 5 x 
10^ time steps to reach a steady state, and the velocity 
profile is then averaged over long time intervals (up to 
10^ time steps) with horizontal bins of thickness Ay = 
0.02(7. Close to the wall, the velocity is zero because no 
particles are present in that region. In the middle of the 
channel, the velocity profile exhibits a parabolic shape, 
such that we can fit the velocity to a plane Poiseuille flow 
and obtain an effective slip length. The fitting function 
we used is 

v{y)=^-^{y\-v' + 2bvB). (53) 

where F*"^' is the external force applied to each DPD 
particle to induce the pressure gradient. For the po- 
sition of the hydrodynamic boundary ys, we use the 
value obtained from simulations for homogeneous sur- 
faces (1.06±0.12(7 from the simulation wall). We assume 
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FIG. 3: Typical density (a) and velocity (b) profiles simulated 
for a longitudinal flow and a texture with L = H = h = 50a, 
(1)2 = 0.5. 

that the position of the hydrodynamic boundary remains 
the same for striped surfaces. The values of {Q)x and 
{Q)z arc then computed with Eqs. ([T3| - (fT4)) . 



V. RESULTS AND DISCUSSION 

In this section, we present the DPD simulation results 
and compare them with predictions from the continuous 
theory. Throughout the section, we use striped patterns 
with periodicity L = 50(t and slip length b = 50(t on the 
slippery areas. We study the effect of varying the film 
thickness H , the angle 8 between the applied force and 
the stripes, and the area fraction <f)2 of slippery areas. 

We start with varying in a system where the stick 
and slip areas are equal, 02 = 0.5, which corresponds 
to maximum transverse flow in a thin channel situation 
[28j . Fig.|3]shows the results for the effective downstream 
slip lengths, , as obtained from Poiseuille fits (|55|) to 
the X component of the velocity Vx, for two values of 
the channel thickness, H = SOtr and H = lOcr. We em- 
phasize that in both cases we formally have an interme- 
diate channel situation, since H/L = 0(1). The error 
bars have been obtained from averaging over five inde- 
pendent runs. Since we model an infinitely extended slit 
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FIG. 4: The effective downstream slip length, h),f( , as a func- 
tion of tilt angle Q for a pattern with L = b — 50a and 
(f)2 ~ 0.5. Symbols are simulation data. Solid lines are theo- 
retical values calculated using Eq. (|54|l with eigenvalues ob- 
tained by a numerical solution of Eqs. (|19|) . (|20|l and (|30|l . 
pip : (a) H — 50a. The thick channel limit (dashed line) is 
calculated with Eqs. ([22} and (b) H = lOcr. The thin 

channel limit (dashed line), is calculated with Eqs. (|2ip and 

by virtue of applying periodic boundary conditions in the 

(x) 

xz plane, the downstream slip length in our system, 6^^^ , 
corresponds to 6°*^. (In channels that are confined in the 
y direction, a transverse pressure builds up that renor- 
malizes 5|,g''[5|.) Fig. |4]also includes a theoretical curve 
calculated using Eq. ([T]), which can be explicitly written 
as 

6f = 6^ffCos2e + 6e'fiSin'e. (54) 

Here, the eigenvalues of the slip-length tensor are ob- 
tained by numerical solution of the dual series, using the 
procedure described in Ref. [2l[ . Also included in Fig. 0] 
(dashed curves) are the xx components of the effective 
slip length tensor calculated in the limit of thick [Fig. 
m^a)] and thin [Fig. IH^b)] channels. These indicate the 
range of b°^^ in this (symmetric) channel geometry. 

The simulation data are in a good agreement with the- 
oretical predictions, confirming the anisotropy of the flow 
and the validity of the concept of a tensorial slip for arbi- 
trary channel thickness [2l|. Most notably, the effective 
slip is larger for thinner channels, somewhat in contrast 
to previous findings for asymmetric channels, where the 
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magnitude of the effective slip length in a thin gap was 
much smaller than that in a thick channels pll l29l [30|. 
This means that transverse hydrodynamic phenomena 
should be enhanced in thinner channels. 
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FIG. 5: The ratio {Q)z/{Q)x as a function of the tilt angle O 
obtained with Eq. ((39| for the data sets from Fig. [4] for (a) 
H = 50a, (b) H = 10a. Symbols are simulation data, solid 
curves represent theoretical values, and dashed curves show 
asymptotic predictions in the limit of thick channels (a) and 
thin channels (b). 

The development of transverse flow is illustrated in 
Fig.[5l where the data for the ratio {Q)z/{Q) x of the mea- 
sured averaged longitudinal and transverse flow rates at 
different are presented and compared with the theoret- 
ical prediction of Eq. (pQ]) . The simulation data and the 
theoretical prediction agree quite well. There are some 
discrepancies at intermediate angles 0, where the sim- 
ulation data for {Q)z/{Q)x are smaller than predicted 
by theory. These slight deviations are most likely a re- 
sult of the finite Reynolds numbers, which are of order 
0(10) in the present simulations as discussed in Section 
IfVI When increasing the bulk force and hence the aver- 
age flow velocities, the deviations increase, thus they will 
presumably vanish in the Stokes limit. A pressure gradi- 
ent in eigendirections cannot produce any transverse flow 
(and the tensorial boundary condition reduces to a scalar 
one), which is well seen in Fig. [5l In all other situations 
the direction of flow is different from that of the pressure 



gradient. The maximum value of {Q)z/{Q)x is smaller for 
the larger channel (Fig.[51[a)). In this case, this maximum 
is reached at 6 = 7r/4, which agrees with theoretical cal- 
culations made with Eq. (^51) for a thick channel. For the 
thinner channel, however, the maximal {Q)z/{Q)x is ob- 
served at larger 8 (Fig.lSfb)), also in agreement with the 
theory (see discussion above). An important conclusion 
from these results is that the surface textures which op- 
timize transverse flow differ significantly from those opti- 
mizing effective (forward) slip. Similar predictions have 
been made for other channel configurations [13, • 

Next we examine the effect of varying the fraction of 
slippery gas/liquid phase, (j)2- Fig. [H] shows the eigenval- 
ues, big and 6^, of the slip- length tensor as a function 
of 02. The two data sets again correspond to different 
thickness of the channel, H = 5Qa and H = 10a. The 
results clearly demonstrate that (f>2 is the main factor de- 
termining the value of effective slip, which significantly 
increases with the fraction of the slippery sectors. The 
theoretical curves match the simulation data very nicely. 




(b) 

FIG. 6: The eigenvalues of the effective slip length tensor 
(symbols) as a function of <j}2. Solid curves are theoretical 
values obtained by a numerical solution of Eqs. (|19p . (|20p 
and (|30p . (|3ip . Calculations were made for a pattern with 
L = b = 50a: (a) H — 50a. Dashed curves are calculated 
with Eqs. ([22]) and ((38l) (b) H = 10a. Dashed curves are 
computed with Eqs. (pT|) and ((35]). 

To illustrate the effect of (/)2 on the transverse phenom- 
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ena, we now fix = 7r/4 and measure {Q)z/{Q)x- The 
simulation data presented in Fig. [7] show that the max- 
imum value of {Q)z/{Q)x is observed at 02 > 0.5 and 
depends on the channel thickness. As discussed above, 
(/)2 = 0.5 results in maximum transverse flow in a thin 
channel (thin compared to the period of the texture) [1^ . 
For thick channels, the continuum theory predicts max- 
imal transverse flow at 02 — 1- The simulation data, 
obtained for two values of H, confirms these trends. Fig. 
[7] also illustrates that transverse phenomena are more ef- 
ficient in a thinner channel. The data sets in Fig. [7] are 
in quantitative agreement with predictions of the contin- 
uum theory. Only for the thick channels and (j)2 > 0.5 
do we observe slight deviations, which again presumably 
reflect inertia effects. 
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FIG. 7: The ratio {Q)z/{Q)x as a function of <f)2 obtained 
with Eq. (|39|) for O = 7r/4 by using the data sets from Fig. 
[6] (a) H = 50(7, (b) H = 10a. Symbols are simulation data, 
solid curves represent theoretical values, and dashed curves 
show asymptotic predictions in the limit of thick channels (a) 
and thin channels (b). 

The simulation and theoretical data presented in Figs. 
U] and |6] show that the eigenvalues of the slip- length ten- 
sor depend on the channel thickness if, in agreement 
with earlier predictions (2l| . To examine this depen- 
dence in more detail, we now fix 02 = 0.5 and vary 
H in a large range. The smallest gap was taken to be 
H = 5ct, for smaller gaps the velocity profile can no 



longer be resolved satisfactorily. The largest thickness 
was chosen to be equal to L = 50a. The eigenvalues of 
the slip-length tensor are shown in Fig. |8l Also included 
are continuum theoretical calculations (solid curves) and 
asymptotic values (dashed curves) expected in the true 
limits of thin and thick channels. The data show that 
the longitudinal effective slip decreases with H. In con- 
trast, the transverse effective slip increases with H, so 
that the difference between two eigenvalues is largest for 
thin channels. We remark and stress that the effective 
slip reaches the asymptotic values predicted for a thick 
channel already at H ~ 0{L). A similar observation was 
made in pTj for asymmetric channels. This result is re- 
markable since it suggests that the thick channel limit, 
where the effective slip is a property of single interfaces 
and does not depend on H, is already reached for chan- 
nels whose thickness is of the order of the texture period. 
In practice, this implies that huge simulations with large 
simulation boxes are not necessary to determine the ef- 
fective hydrodynamic behavior in such systems. We note 
however that even for the smallest gaps considered here 
{H/L = 0.1), our theoretical and simulation results still 
deviate strongly from the upper Wiener bound, Eq. (PT|) . 
although they are very close to the lower Wiener bound, 
Eq. dSSD. 

Fig. ini shows the corresponding simulation and theo- 
retical data for {Q)z/{Q)x as a function of channel thick- 
ness H. The stripe texture is the same as in Fig. [51 
i.e., 02 = 0.5, and the tilt angle was chosen 9 = n/A. 
As discussed above, this angle leads to a maximal trans- 
verse flow in case of a thick channel, but not in the thin 
channel situation, where a much larger angle is required 
to optimize the transverse phenomena. Nevertheless, 
{Q)z/{Q)x increases dramatically, when the channel be- 
comes thinner. With our range of parameters, however, 
we are still well below the limiting value of {Q)z/ {Q)x, 
which would be expected in case Wiener bounds were 
attained. 

Finally in this section, we show the measured veloc- 
ity profile in a channel with stripes inclined at the angle 
= 7r/4 with respect to the applied force fFig. [TO)) . The 
velocity profile results from superimposed fiows in the x 
(forward) and the z (transverse) directions (also shown). 
As above we use H = lOtr [Fig. [TUTa)] and H = 50a 
[Fig. [TOlb)]. The data show that the effective velocity 
direction is generally misaligned with the force vector, 
and this effect is much more pronounced for a thin chan- 
nel due to a larger transverse component of the velocity. 
The flow has the form of a superimposed no-slip parabolic 
Poiseuille flow and a slip-driven plug-fiow everywhere in 
the thin channel. In case of a thick channel the fiow 
near the surface is different from that in the center of the 
channel. To examine this more closely, the short-distance 
region of the flow from Fig. [TUT b) is reproduced in Fig. 
[TUT c). and the central part is shown in Fig. [TUfd). These 
blowups demonstrate that the strong transverse flow due 
to surface anisotropy is generated only in the vicinity of 
the wall and tends to disappear far from it, as predicted 
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FIG. 8: The longitudinal (a) and transverse (b) effective slip 
lengths as a function of the channel height H for a texture 
with L = b = 50a and <f)2 = 0.5. The theoretical curves are 
obtained by a numerical solution of Eqs. H19|) . (|20|) and (|30|) . 
(|31|) . Dashed lines show expected asymptotics in the limit of 
thin and thick channels. 
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FIG. 9: The ratio between the transverse and longitudinal 
flow rates {Q)z/{Q)x as a function of channel thickness H for 
a pattern with L — b = 50a, 02 = 0.5, and — n/4. Sym- 
bols are simulation data and the lines represent theoretical 
values obtained using Eq. p9|l . Dashed lines show expected 
asymptotics in the limit of thin and thick channels. 



theoretically [T3l and observed in experiment H^]. The 
effective velocity profile is 'twisted' close to the striped 
wall. 



VI. CONCLUSION 



We have investigated pressure-driven flow in a flat- 
parallel channel with two aligned striped SH surfaces. 
For this geometry of configuration wc propose a semi- 
analytical theory, which is valid for an arbitrary gap and 
any local slip length at the slipping areas of the stripes. 
Analytical results are presented for various important 
limits, and our approach gives simple analytical formula 
for an effective slip length in case of stripes that are in- 
clined with respect to a pressure gradient. Furthermore, 
our theory gives analytical guidance as to how to choose 
the parameters of the texture and a tilt angle, in order to 
optimize the transverse flow in different situations. Our 
theoretical predictions have been compared with results 
of DPD simulations, which are in a good quantitative 
agreement with the macroscopic theory. Our results are 
directly relevant for passive mixing in thin channels and 
other microfluidic applications. In the future, we hope 
to augment our analysis to include the possibility of a 
surface charge patterns, which would help to understand 
transverse electro-osmotic phenomena. Another fruitful 
direction could be to consider identical, but misaligned 
striped walls, where the relative orientation of surfaces 
could generate additional mechanism of mixing. 
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